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Prediction  of  Materials 
Charging  in  Magnetospheric  Plasmas 


I.  INTRODUCTION 

It  is  often  useful  to  obtain  a  quick  approximation  of  the  potential  to  which  a 
material  would  charge  in  a  geosynchronous  orbit  plasma.  If  one  does  not  con¬ 
sider  multidimensional  effects  and  the  effects  of  photoelectrons,  a  computation 
of  the  materials  effects  alone  is  possible. 

This  calculation  is  carried  out  by  means  of  a  computer  code  called  MATCHG 
(Materials  Charging).  Given  a  plasma  environment  of  specified  electron  and 
proton  temperature  and  number  density,  the  code  calculates  the  potential  to  which 
a  material  will  charge  in  this  plasma.  The  principle  of  the  method  is  to  solve  for 
the  potential  for  which  there  is  zero  net  current  to  the  material  surface.  Included 
in  the  current  balance  are  incident  electron  and  proton  (ion)  currents,  secondary 
electrons  from  both  incident  electrons  and  ions,  backscattering,  and  current 
leakage  to  a  substrate  through  bulk  conductivity.  In  the  formulation  which  includes 
bulk  conductivity,  the  thickness  of  the  material  must  be  specified.  The  secondary- 
emission  yield  curves,  as  a  function  of  energy,  are  calculated  in  the  code. 


(Received  for  publication  17  March  1981) 


5 


2.  MATERIALS 


The  materials  for  which  calculations  were  carried  out  are  those  typically 
used  on  satellite  surfaces  (see  Table  1).  "SOLAR"  refers  to  Solar  cell  cover 
glass  material;  "WHITEN"  to  a  nonconducting  white  pah.i;  "  \  Hi  <  )WC"  to  a  1! 
cotiiltie’.ng  paint;  "BLACKC"  to  a  black  conduct  paint. 

1  he  MATt'llG  pt  ,  gran:  barges  a  rftiiT; .  by  mean-;  .  f  tr.aterial  "  "- 

•t.uir.t  ;•»;>&  :  ii.it  are  used  in  NASC’A  P.  *  ;(  •  a  -'  1,  :  one  •  o’  tt.irU  •  r 

the  equilibrium  potentials  for  a  surface  immersed  in  a  spe  <1  plasma  envir:  n 
ment.  The  full  3-D  complexities  of  more  realistic  models  are  neglected. 

In  this  report,  we  present  the  main  results  from  MATCHG.  Graphical  and 

linearized  techniques  are  developed  for  approximating  the  results  of  the  MATCHG 

equations.  Runs  were  made  with  the  plasma  temperature  T  (T.  -  2T  )  as  the 

e  t  e 

independent  variable.  The  densities  were  fixed  at  np  =  n.  =  10  M  .  The  poten¬ 
tials  were  computed  in  three  ranges; 

1.  Low  temperatures  (see  Figure  1) 

2.  Intermediate  temperatures  (see  Figure  2) 

3.  High  temperatures  (see  Figure  3). 

Runs  were  made  to  calculate  the  potential  for  a  two-Maxwellian  plasma  case 
and  the  corresponding  single  Maxwellian  plasma  case.  The  parameters  specify 
three  types  of  plasmas:  hot,  cold,  and  average.  These  are  shown  in  Table  2. 

The  MATCHG  results  are  given  in  Table  3  for  Teflon  and  aluminum. 

The  intermediate  temperature  cases  were  done  for  all  of  the  materials  cur¬ 
rently  available  in  MATCHG  (see  Figure  2).  Note  that  materials  BLACKC  and 
YELOWC  are  identical,  that  is,  they  are  represented  in  the  program  by  the  same 
material  parameters.  Also  the  materials  WHITEN  and  Kapton  are  sufficiently 
near  to  YELOWC  and  BLACKC  that  the  final  potentials  are  the  same  (to  three 
decimal  places). 


3.  CHARGING  AT  COW  TEMPERATURES 

Spacecraft  charging  is  most  significant  for  plasma  temperatures  greater  than 
1  kv.  It  is  of  interest  to  study  charging  at  lower  plasma  temperatures  as  well. 

For  nonsunlit  conditions,  Figure  1  shows  the  equilibrium  potentials  for  a 
number  of  materials,  both  metals  and  dielectrics. 


1.  Katz,  I.  et  al  (1977)  A  Three  Dimensional  Dynamic  Study  of  Electrical  Charg¬ 
ing  in  Materials.  NA&A  CR-135256. 
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Table  2.  Double  and  Single  Maxwellian  Plasma  Cases 


l) 

Hot  Plasma 

n  ,  =  0.  06 
el 

n  ,  =  0.06 
i  1 

"e2  =  0. 14 

"i  2  '  °‘14 

T  ,  -  1200  eV 
el 

T  -  2000  eV 
il 

T  „  =  10000  eV 
e2 

T...  =  16000  eV 

i2 

1 -Maxwellian  version: 

n  =  n.  =  0.  20 
e  i 

T  =  8000 
e 

T.  =  13000 

i 

2) 

Cold  Plasma 

n  ,  =0.85 
el 

n.,  =  0.25 
il 

ne2  =  0.25 

"i2  =  °'85 

T  ,  =  1000  eV 
el 

T.,  =  300  eV 
i  1 

T  9  =  4000  eV 

T. _  =  9000  eV 

i2 

1 -Maxwellian  versions: 

n  =  n.  =  1.  1 
e  i 

T  =  1800 
e 

T.  =  8000 

i 

3) 

Average  Plasma 

"el  =°-9 

"ii  =°-9 

"e2  =°-6 

"i2  °-  6 

T  ,  =  1500  eV 
el 

T. ,  =  427  eV 
il 

Tp2  =  4400  eV 

T(2  =  10000  eV 

1 -Maxwellian  versions: 

n  =  n.  =  1.  5 
e  i 

T  =  2800 
e 

T.  =  7000 

i 

V  (VOLTS) 
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Figure  1.  AlATCHG  Potentials  for  Low  Temperatures 


T,  <k  «V) 


Figure  2.  MATC'HG  Potentials  for  Intermediate  Temperatures 


Figure  3.  IUATOHG  Potentials  for  High  Temperatures 


At  extremely  low  temperatures,  up  to  5  eV,  all  materials  charge  negatively 
to  approximately  the  same  potential.  Between  10  and  100  eV,  secondary  emis¬ 
sion  properties  become  important  and  the  potentials  of  different  materials  peak 
between  -20  V  for  Teflon  to  -115  V  for  gold.  Between  about  20  eV  and  200  eV 
plasma  temperature,  different  materials  charge  to  positive  voltages.  The  posi¬ 
tive  voltages  are  all  less  than  5  eV,  most  around  2  eV. 

In  the  temperature  region  between  10  and  several  hundred  eV,  materials 
potentials  are  either  a  few  volts  positive  or  tens  of  volts  negative. 

These  results  are  significant  in  the  design  of  electron  spectrometers  where 
the  potentials  produced  by  low-energy  electrons  impinging  on  dieleetric  materials 
is  important. 

4.  FFFF.OT  OF  BM.K  COMM CTIVITY 

IUATCIIG  includes  the  effect  of  the  hulk  conductivity  current.  In  order  to 
calculate  this  current,  the  potential  difference  across  the  charging  surface  must 
be  known,  that  is,  the  voltage  underlying  the  surface  is  needed. 
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Table  3.  Single  vs  Double  Maxwellian  Results 


In;t  tal 

Final 
Potent  tal 

K  :it  i  i  >.s 
(MAX  ' 

Kle.  W‘.,M 
(  urn-nt 

l:  •  :  ... 

(MAN 

Materia! 

(keV) 

1  till  .MAX  1 

t!0  , r:  -1 

DPI  ..M  Ac.  i 

MAN  -  h  ' 

Aluminum 

-r>.  oo 

1.01 

-o.  47'*  •) 

’  1  • 

DUI  MAN  -  !.  > 

Aim: 

-M  M  J 

-o4  ;  o J 

MAN  -  a\  .■ 

Alumunm 

-1.11  1 

in  - 

0.  «« 

'} 

Dill  .MAN  -  ave 

Aluminum 

-1.-.4  J 

-2. 00  J 

MAN  -  void 

Aluminum 

-o.i I 

l.  00 

-m,  , 

1  .  1(1 

DIM. .MAN  -  cold 

Alum  mum 

-0.  oho r,  1 

-1.  14  > 

MAN  -  hot 

Teflon 

-,?2  1 

0.  04 

-0.470  I 

1.11 

DRI.MAN  -  hot 

Teflon 

— v 

o 

ei 

i 

-0. 420  > 

MAN  -  ave 

Teflon 

-0. 00071  l 

l.  0«; 

1.00 

DRI.MAN  -  ave 

Teflon 

0. 00007  / 

-2.00  > 

MAN  -  t  old 

Teflon 

+0.  00127  I 

0.  92 

1.  10 

Teflon 

+0.  00128  / 

-1.14  J 

To  make  a  comparison  to  a  NASCAP  satellite  model,  the  potential  on  the 
underlying  conductor  would  be  taken  as  the  backing  plate  voltage.  Table  4  shows 
the  shade-side  potentials  for  quasispheres  in  a  fixed  sun  environment  and  the 
corresponding  MATCHG  voltages,  including  the  bulk  conductivity  current.  There 
is  good  agreement  for  all  of  the  materials. 

In  Figure  4  we  plot  the  MATCHG  potentials  vs  the  backing  plate  voltage  for 
the  case  of  Teflon  in  the  nominal  T  =  10,  20  keV  plasma.  This  figure  shows  that 
the  minimum  voltage  |  -v|  occurs  when  the  backing  voltage  is  0. 

In  Figure  5  we  show  the  lower  limit  MATCHG  potentials  (backing  plate  voltage 
set  to  0)  for  a  number  of  different  materials.  The  plasma  temperature, 

T  (T.  =  2T  ),  is  the  independent  variable  and  the  densities  are  fixed 
(np  =  n.  =  10  M  ).  It  is  evident  that  those  materials  with  large  bulk  conductiv¬ 
ities  (WHITEN,  SI02)  are  depressed  relative  to  the  others,  and  SOLAR,  with  a 
small  conductivity,  charges  up  higher,  that  is,  the  larger  the  conductivity,  the 
easier  it  is  for  charges  to  leak  through  the  material  and  reduce  the  voltages. 


11 


Table  4.  MATCHG  Potentials  vs  NASCAP  Quasisphere 


Material 

NASCAP 

Shade 

Potential 

Quasisphere 

Conductor 

Potential 

MATCHG 

Backing 

Plate 

Final 
Potent  ial 

Teflon 

-8894 

-7125 

-7100 

-8800 

Kapton 

-9210 

-7378 

-7380 

-9280 

SOLAR 

-7820 

-02  04 

-0200 

-7820 

SI02 

-358 

-279 

-280 

-300 

WHITEN 

-2997 

-2395 

-2390 

-2990 

Plasma:  T  -  10,  T  =  20  keV,  n  =  n.  =  10b  M 
e  '  i  'ei 

Sun  is  in  direction  (1,  1,  1) 


Figure  4.  MATCHG  Potentials  vs  Backing  Voltage  (Bulk  Conductivity 
Included) 
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Figure  5.  MATCHG  Potentials  (Bulk  Conductivity 
Included) 


Lower  limit  curves  of  this  type  would  be  directly  applicable  to  cases  where 
the  spacecraft  conductor  does  not  charge  and  the  surface  lies  in  the  shade  (and 
the  surface  conductivity  is  negligible).  For  example,  we  can  compare  the  num¬ 
bers  obtained  for  the  SCATHA  model  for  experiment  SC  1-3  (which  lies  in  the 
shade)  to  the  MATCHG  results,  as  shown  in  Table  5.  The  agreement  is  very 
good. 

When  the  bulk  conductivity  current  is  included  in  the  computations,  the  steady- 
state  potentials  are  determined  by  a  balance  between  the  current  arriving  from 
the  plasma  and  that  escaping  through  the  surface  to  the  backing  plate.  The  final 
voltages  reached  will  thus  depend  on  the  plasma  density  (which  is  related  simply 
to  the  plasma  current).  This  was  not  the  case  before  when  the  bulk  conductivity 
was  neglected  (see  Figures  1  and  2),  In  Figure  6,  we  show  the  large  changes 
that  take  place  in  the  MATCHG  potentials  when  the  plasma  density  is  varied  by 
factors  of  10  and  0.  1  (for  the  case  of  a  Teflon  surface). 

One  can  get  a  feeling  for  the  effect  of  the  ion  temperature  of  the  plasma  on 
the  MATCHG  potentials  from  Figure  6.  The  plot  shows  the  change  in  the  voltage 
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curve  that  occurs  when  we  use  the  relation  Tg  =  T.  instead  of  the  nominal 
Te  =  2T;.  There  is  only  a  small  increase  (~  10  percent)  in  the  equilibrium  poten¬ 
tials. 

An  alternate  formulation  for  the  secondary  emission  was  added,  as  a  test 

mode,  to  MATCHG.  In  this  new  version,  the  Sternglass  yield  equation  is  used 
2 

at  low  energies 


By  imposing  continuity  of  slope  and  value  for  the  above  two  equations  we  set 
the  transition  energy  x  =  3  and  the  ratio  K/im  =  1.  54.  The  MATCHG  potentials, 
from  runs  using  this  test  mode,  are  shown  in  Figure  6.  The  test  equations  result 
in  a  slight  increase  in  the  number  of  secondaries  and  thus  the  equilibrium  poten¬ 
tials  are  decreased,  as  indicated  in  the  figure. 


5.  GRAPHICAL  AND  LINEARIZED  REPRESENTATION  OF 
CHANGING  EQUATIONS 

We  solve  the  circuit  equations  representing  the  charging  calculation  under 
the  assumption  that  the  source  of  current  is  a  linear  function  of  the  source  volt¬ 
ages.  This  procedure  gives  a  fairly  good  approximation  to  the  MATCHG  results  — 
the  general  features  of  the  charging  process  are  obtained  (including  final  equilib¬ 
rium  voltage  and  time  constant,  which  depends  on  conductivity,  thickness  of 
material,  etc. ). 


6.  CIRCUIT  EQUATIONS 

The  charging  calculations  that  are  done  in  the  MATCHG  program  can  be  rep¬ 
resented  by  the  circuit  shown  in  Figure  7.  The  equation  describing  the  voltage 
on  the  capacitor  is,  from  the  diagram, 

2.  Dekker,  A.J.  (1958)  Solid  State  Physics,  Academic  Press,  New  York  8:279. 

3.  Wall,  J.A.  et  al  (1977)  Results  of  Literature  Search  on  Dielectric  Properties 

and  Electron  Interaction  Phenomena  Related  to  Spacecraft  Charging, 
Proceedings  of  the  Spacecraft  Charging  Technology  Conference.  AFG1.- 
TH  -77-0051,  AD  A045  459. 
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dV 

dt 


i  -  (V  -  Vfi)a  A/d 


Ke  A/d 

o 


1 

Ke, 


(dj  -  (V  -  VB)a) 


(1) 


where  j  i/A  is  the  new  current  density  from  the  source  (it  does  not  include  the 
bulk  conductivity  current).  MATC'HG  calculates  j,  which  is  a  complicated  func¬ 
tion  of  the  plasma  environment  and  the  properties  of  the  surface  material.  In 
Figure  8,  we  show  the  j  j(V)  dependence  for  Teflon  in  a  plasma  specified  by 

T  5  keV  T.  =  G,  n  0.  1  cm  ,  n.  =  0.  15.  It  is  evident  that  the  i  vs  V  curve 
e  ■  ■  i  ’  e  t  J 

is  almost  linear.  Thus  it  will  be  a  fairly  good  approximation  to  take 


j  =  j 


+ 

o 


J'V  , 


where  J'  is  some  average  slope.  For  simplicity,  we  could  let 


(2) 


J' 


V(o  =  0)  ’ 


(3) 


where  j  =  initial  net  current  density,  V(a  =  0)  =  equilibrium  potential  when  a  =  0 
which  gives  the  straight  line  indicated  in  Figure  8  (V(o  =  0)  =  0  and  j  <  0  so  that 
J'  <  0).  Substituting  Eq.  (2)  into  Eq.  (1)  gives 

dF=  i<7  (dj0  +  vB)  +  v(iA-(JI -°/d))=  A  +  BV  •  (4) 

O  O 

where  we  define 

A  =  K7-<d)0+VBa) 
o 

B  =  -  a/d)  <  0  . 

O 

Equation  (4)  is  easily  solved.  We  write  it  as 

(a'+BV)  =  dt  •  (7> 

Integrating  both  sides  of  Eq.  (7)  gives 

^ln  (A  +  BV)  =  t  +  C  .  (8) 


(5) 

(6) 
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*•  ^dielectric  constant  dv^jc 

<0*8.85  *  I0'12  dt  *  C 

d  *  thickness  of  material  ,  . ,  _ , 

A = AREA  OF  SURFACE 
a-  *BULK  CONDUCTIVITY 


Figure  7.  MATCHG  Circuit  Diagram 


EQUATIONS 
FOR  CAPACITOR 
CHARGING 
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The  integration  constant  C  can  be  evaluated  by  letting  t  =  0  in  Eq.  (8),  that  is, 


C  =  4- ln  (A  +  B  V  ) 

r$  O 


(9) 


where  Vq  =  V(t  =  0).  Putting  Eq.  (9)  into  Eq.  (8)  and  solving  for  V  leads  to 


(10) 


In  Figure  9  we  plot  V  from  Eq.  (10)  (for  the  case  Vq  =  0)  vs  the  time  and  the 
corresponding  results  obtained  from  MATCHG.  The  predicted  curve  is  slightly 
higher.  The  time  constant  for  charging  up  can  be  taken  as 


1 

(J'  -  o/d) 


Keo 


1 

(J'd/a  -  1) 


(11) 


Check:  when  J<  =  0  (constant  source  current),  we  have 


which  is  what  we  would  expect  from  elementary  circuit  theory.  The  final 
equilibrium  potentials  can  be  obtained  by  setting  =  0  in  Eq.  (4),  or  from 
Eq.  (10)  by  setting  t  =  oo,  since  the  exponential  term  goes  to  0  (B  <  0). 


Ve  =  V(«o)  = 


A_ 

B 


-<jp  +  VB  °/d> 
(J'  -  a/d) 


or  using  Eq.  (3), 


=  V(ct  =  0) 


Check:  when  a  =  0  we  recover  the  result  that  V  =  V(ct  =  0). 

e 


(12) 
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-  -657  =  CALCULATED  LIMIT - 


-609'MATCHG  LIMIT 


AfeB,-l) 
B  \ 


'  MATCHG 


0  1000  2000  3000  4000  5000 

t  (sec)—*' 

Figure  9.  Calculated  V  vs  Time  From  MATCHG 


7.  GRAPHICAL  SOLUTION 


As  time  increases,  the  source  current,  i,  decreases  and  the  conductivity 
current,  i^,  increases  in  proportion  to  the  voltage  buildup  on  the  capacitor.  At 
equilibrium,  the  current  to  the  capacitor,  i^-.,  has  decreased  to  0  and  thus 
i  =  ijj,  or  in  terms  of  current  densities  we  have 


i  in  +  J'Vo  =  j| 


—  (V  -  v  ) 
d  e  VB' 


Solving  the  condition  for  Vg  again  yields  the  Eq.  (12).  One  could  solve  the 
above  equations  graphically  by  finding  the  intersection  of  the  straight  lines  rep¬ 
resenting  j  and  jp  on  the  (j  ,  V)  plot  (see  Figure  8).  In  fact,  an  exact  solution 
would  be  determined  by  the  intersection  of  the  ,i^  line  and  the  MATCHG  j(V)  curve. 
The  effect  of  changing  a,  d,  and  Vg  is  clear  from  this  construction.  One  just 
changes  the  slope  and  intercept  of  the  line  and  finds  the  new  intersection  point 
giving  Ve. 
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8.  NUMBERS  FORT  =  5.6  CASE 


Kapton 


-14 


a  -  10 

Via  --  0)  =  -4.43  X  10 


K  =  3.  5 
3 


(-5.  83  X  108) 
(-4.43  X  10'3) 


d  =  1.  27  X  10 


-4 


i  =  -a.  83  X  10 
Jo 


-8 


=  -1.3  X  10 


-11 


A  = 


B 


dj, . 


(1.  27  X  10‘4)(-a.  83  X  10"8) 


Ke 


(3.  5)(8.  85  X  10'12) 


=  k7~  <Jt  *  CT/d> 


1.  27  X  10 


-4 


(3.  5)(8.  85  X  10‘12) 


-1.3  X  10 


-11 


=  3.  64  X  10 


-4 


-0.  239 


10 


-14 


1.  27  X  10 


-A 

B 


(-0.  239) 


(-3.  64  X  104) 


=  -657 


Ke, 


(*•1-0 


(3.  5)(8.  85  X  10'12) 


10 


rnr 


(-1.3  X  10'n)(1.27  X  10 ~4 


10 


r 


2.  66  X  10J 


=  -609. 


from  MATCHG  run,  V 


IATE 

LMED 


